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ABSTRACT 

We explore the properties of cold gravitoturbulent accretion disks — non-fragmenting disks hovering on 
the verge of gravitational instability — using a realistic prescription for the effective viscosity caused 
by gravitational torques. This prescription is based on a direct relationship between the angular 
momentum transport in a thin accretion disk and the disk cooling in a steady state. Assuming that 
opacity is dominated by dust we are able to self-consistently derive disk properties for a given M 
assuming marginal gravitational stability. We also allow external irradiation of the disk and account 
for a non-zero background viscosity which can be due to the MRI. Spatial transitions between different 
co-existing disk states (e.g. between irradiated and self-luminous or between gravitoturbulent and 
viscous) are described and the location of the boundary at which disk must fragment is determined 
in a variety of situations. We demonstrate in particular that at low enough M external irradiation 
stabilizes gravitoturbulent disk against fragmentation all the way to infinity thus providing means of 
steady mass transport to the central object. Implications of our results for the possibility of planet 
formation by gravitational instability in protoplanetary disks and star formation in the Galactic Center 
and for the problem of feeding supermassive black holes in galactic nuclei are discussed. 

Subject headings: accretion, accretion disks — instabilities — (stars:) planetary systems: protoplan- 
etary disks — (galaxies:) quasars: general 



1. INTRODUCTION. 

Gravitational instability (GI) in astrophysical disks has 
been a subject of investigation for more than fifty years 
since the seminal work by Safronov (1960) and Toomre 
(1964). Originally it was studied predominantly in the 
context of driving the spiral structure in galaxies. Later 
it has been suggested that GI may play important role 
in planet formation (Cameron 1978; Boss 1998), and its 
significance for the properties of compact nuclear disks 
around supermassive black holes in the centers of galaxies 
has also been recognized (Paczynski 1978a; Goodman 
2003). 

It is generally accepted that the GI sets in when the 
so called Toomre Q satisfies the following condition: 



Q 



fit 



tuations but the disk is stable against fragmentation in 
the long run and maintains itself on the brink of instabil- 
ity with Q w Qq. Stability against fragmentation arises 
because the restoring action of the thermal pressure re- 
sisting the self-gravity of overdensities is not sufficiently 
reduced by cooling when t coo i > O . 

Torques produced by the nonaxisymmetric density per- 
turbations in gravitoturbulent disks give rise to angular 
momentum transport. Considering the disk to be in a 
steady state on time shorter than its viscous evolution 
timescale one can directly relate "effective viscosity" aci 
produced by the gravitational torques to the cooling time 
of the disk. Gammie (2001) has demonstrated that the 
dimensionless a-parameter characterizing disk viscosity 
(Shakura & Sunyaev 1973) is 
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where £, f2, and c s axe the local surface density, angu- 
lar frequency, and the sound speed in the disk which we 
consider to be made of gas and having a Keplerian ro- 
tational profile. Qo is a constant of order unity, its pre- 
cise value determining the instability threshold ranges 
from 0.7 to 1.7 according to different authors (Kim et 
al. 2002; Boss 2002). Nonlinear development of the GI 
sensitively depends on the thermodynamical properties 
of the gas as has been first shown by Gammie (2001): if 
gas can cool on a timescale shorter than the local dynam- 
ical timescale ft^ 1 then disk fragments into bound, self- 
gravitating objects which may grow further by accreting 
the surrounding gas and colliding with each other. But 
if the cooling timescale t coo i is longer than f2 _1 then the 
disk settles into a state of the so-called gravitoturbulence 
in which surface density can experience significant flue- 
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It is instructive to show where this relation comes from. 
First, in a steady state the rate of energy dissipation per 
unit surface area of the disk ~ £1 2 M has to equal the 
energy flux F emitted from the disk surface. Second, the 
accretion rate is M = 37rz/E, where v = ac^/ft (Pringle 
1981). Combining these relations one immediately ob- 
tains equation ^ with 



cool 



Ec 2 s /F. 



(3) 



Relation ([2|) makes it possible to interpret disk fragmen- 
tation occurring at t CO oi ~ Si as the inability of the 
disk to sustain gravitational stress at a > 1 (Rice et al. 
2005). 

In this paper we investigate the structure and evolution 
of gravitoturbulent disks in which angular momentum is 
transferred predominantly by the gravitational torques. 
This problem has been previously investigated by Lin 
& Pringle (1987) but with a rather naive prescription 
for the effective viscosity. Also, some efforts have been 
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devoted to understanding the structure of the gravito- 
turbulcnt disks which are unstable to fragmentation on 
large scale, i.e. disks having Q ss Qq and otci ~ 1 ev- 
erywhere (Rafikov 2005, 2007; Matzner & Levin 2005). 
In this work viscous evolution of the disk is explored ac- 
cording to the prescription ([2} without fixing the value of 
aai — instead it is calculated self-consistently based on 
the physical properties of the gas. We concentrate our 
attention on rather cool disks in which opacity is due to 
dust grains thus focussing on the GI in the outer parts of 
protostellar disks and disks around supermassive black 
holes. 

2. GENERAL CONSIDERATIONS. 

We consider a gravitoturbulent disk in which the dis- 
sipation of transient density waves excited by GI is ca- 
pable of maintaining Q = Qq, and the cooling time t coo i 
is longer than f2 _1 . Cooling time of the disk is 
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where E is the surface density of the disk, c s = {k B T/ p) 
is the isothermal sound speed determined by the mid- 
plane temperature T, and /(r) is a function of the op- 
tical depth t = J npdz (n and p are the gas opacity 
and density, z is the vertical coordinate) which links the 
emitted flux F to T: F = aT 4 / /(r). A specific form of 
/(r) depends on the way in which energy is transported 
from the midplane of an optically thick disk to its pho- 
tosphere where it is radiated to space. Rafikov (2007) 
has calculated /(r) in the case of efficiently convecting 
disks. However, in this work we assume (as was previ- 
ously done in Rafikov 2005) for simplicity that energy 
is carried from the disk midplane to its surface solely 
by radiation in which case /(r) can be reasonably well 
approximated by 

/(r)«r+i (5) 

T 

This expression smoothly interpolates between the cool- 
ing rates applicable in the optically thick (t > 1) and 
optically thin (r <C 1) regimes. 

We assume a temperature-dependent opacity in the 
form 

k = kqT 13 , (6) 
which is appropriate at low temperatures when n is domi- 
nated by dust grains. At very low temperatures, T < 150 
K, it is generally found (Bell & Lin 1994; Semenov et al. 
2003) that opacity is due to the icy grains and is charac- 
terized by 

(3 = 2 and n a w 5 x l(T 4 cm 2 g -1 K- 2 (7) 
within a factor of 2 or so. At higher temperatures ices 
evaporate and opacity behavior can be crudely described 
as k » 0.1T 1 / 2 cm 2 g" 1 (Bell & Lin 1994). For simplicity 
in this work we do not distinguish between the Rosseland 
mean and the Planck mean opacities (appropriate for 
r> 1 and t«1 correspondingly) as they have similar 
values at low T. 

Now, using definition ([1]) and condition Q = Qo we 
find that 
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in a gravitoturbulent disk. In the optically thick regime 
total optical depth is dominated by the midplane layers 
of the disk in which most of the mass is concentrated, 
so that up to factors of order unity r ~ k(T)E. Clearly, 
this approximation also works in the optically thin case. 
Thus, using equation ([9]) one rather generally finds that 
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We can also calculate aci characterizing angular mo- 
mentum transport caused by the non-axisymmetric sur- 
face density perturbations. Using equations ([2]), ([JJ, and 
J9]) one finds that 
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(11) 



The kinematic viscosity v 
the following expression: 



c*G/c 2 /f2 is then given by 



vgi = C 
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Parameter Q ~ 1 appearing in these equations absorbs 
our ignorance of the exact values of constant factors in 
equations (J2J> and (|4]). 

Equations (fT0|) . (fTT|) and (fl2|) provide us with the de- 
sired viscosity prescription needed for determining the 
physical structure and evolution of the gravitoturbulent 
disk. The only two essential ingredients that went into 
deriving this viscosity recipe are (1) requirement that 
disk maintains itself in a state of marginal stability with 
respect to GI and (2) prescription {2J for aci which 
arises from reasonable assumption that disk is in ther- 
mal equilibrium on timescales shorter than the viscous 
timescale. 

3. CONSTANT m DISKS. 

In this section we consider the structure of gravitotur- 
bulent disk with a specified mass accretion rate M. In 
a steady state M — 37wE, which with the aid of equa- 
tion (Ti"2"l) can be manipulated into the following general 
relation: 



M = 3ttC 
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Function /(t) entering this expression depends on E and 
fl through equation (| 10[) . This allows us to uniquely 
express E as a function of f2 for a given M . 

We should note here that although in the following we 
will mainly discuss disks with constant M our results are 
also directly applicable to disks in which M varies with 
distance. Indeed, as long as one knows M at a partic- 
ular distance (or Q) in the disk equation Q13p uniquely 
determines the value of E in this location 3 . Thus, all our 
subsequent numerical estimates would apply also to the 
case of non-constant M disks as long as M is specified 
at a location of interest. 

3 Strictly speaking the relation M ~ uT, used in deriving (1 1 3 I t 
is valid only in disks with smoothly varying E (e.g. in disks with 
power law dependence of M on r); constant factor in this relation 
is in general different from 37r. 
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We separately consider the cases of optically thick and 
optically thin gravitoturbulent disk. Before we do this 
we note that there are two important transitions charac- 
terizing such disk. One is the 

t = 1 (14) 

transition between the optically thick and optically thin 
regions. Another is the point at which t coo i becomes 
comparable to f2 _1 and disk fragments. This transition 
is defined by condition 



OiGI = X ~ 1, 



(15) 



where x 1S the constant of order unity, its precise value 
has been determined by Gammie (2001) and Rice et al. 
(2003, 2005) in a variety of circumstances. 

According to equations (fTO]) and (fTT|) each of these 
two relations sets a unique constraint between E and il. 
However, if we demand that both of them are fulfilled 
simultaneously (i.e. disk fragments exactly at the r = 1 
transition) then these relations hold only for specific val- 
ues of E and Q, which we denote E/ and Qf. Equations 
(fTO]) . (fTTT) . (fU]) . and (15]) then yield the following values 
of these parameters: 
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(we set /(t) ~ 1 at r = 1). Through equations ©, ©, 
(QjJ) Yjj and Qf (with /(r) ~ 1) also determine fiducial 

values of the mass accretion rate Mf, sound speed c s j, 
and midplane temperature Tf. 
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Note that all these fiducial quantities depend only on 
physical constants and opacity parametrization. For k 
given by equation ([7]) we find the following values of these 
parameters (assuming Qo « 1, Xj C ~ 1) : 



E/ « 15 g cm 
% w 1.4 x 10" 



10 1 
s , 



0.22kms _i , T f 
Mf « 7 x 10~ 6 M 



11.6 K, 
(21) 



The numerical value of Tf conveniently falls into the 
regime of k dominated by icy dust grains. 

For a given mass of a central object M+ angular fre- 
quency f2/ determines a fiducial distance r/ according to 
the formula 
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1/3 



130 AU 



M* 



1/3 



(22) 



In the case of a disk around M+ = 10 6 M & black hole 
one finds Tf ~ 0.06 pc. 

We now consider disk structure in the constant M case 
separately for the optically thick and the optically thin 
regime, as well as the effects of external irradiation and 
the background viscosity in the disk. 

3.1. Optically thick case. 

In the optically thick case /(r) rj t. Plugging this 
into equation (fl3|) and using expression (TIT)]) we find the 
following scalings: 

S = E / (m^ 10 - 2 ^) (7_2,3r 1 = E/m 1 / 3 ^ 2 , (23) 
T = Tf (m 2 ^) {7 - 2Pyl = T/m 2 /^ 2 , (24) 

where second equalities are for (3 = 2 and we have defined 
the following dimensionless quantities: 



M 

m = — — , lu 
M f 



Using equations (p~0|) and (1 1 1 1) we also find that in the 
optically thick regime 

a G z = x(m 4 - 2 ^- 9 ) (7 - 2/3)_1 =x-" 3 , (25) 

t- = ( r | 1 2/3+l,„10+4/3^(7-2/3)- 1 = ^5/3,„6 



(26) 



According to equation (|26[) our assumption of r > 1 is 
self-consistent only if 



uj > u-i = m ( 2/2H 



l)/(10+4/3) 



(27) 

(we assume (3 < 7/2 as expected for dust opacity). Thus, 
gravitoturbulent disk is going to be optically thick for all 
Q > flf (or u > 1) if M > Mf (or m > 1). 

Note that equations ([23)) - (26)) predict very rapid vari- 
ation of disk properties with radius in the optically thick 
regime. Indeed, for /3 = 2 one finds E,T oc r~ 3 , while 
r oc r~ 9 . It is clear that even a moderate increase in 
r would lead to the disk becoming optically thin. Also, 
aci oc r 9 / 2 (analogous scaling has been found by Good- 
man [2003] under different assumptions) and given this 
rapid variation it is obvious that an optically thick grav- 
itoturbulent region can exist only within a limited range 
of radii, see §gfS3& [3331 

3.2. Optically thin case. 

In the optically thin regime we need to use f(r) r _1 
in equation 1T3]) which combined with equation (fTQ|) gives 
the following scalings: 



(28) 
(29) 



T = T f (mc) 2/(9+2/3) = T f m 2 ^co 2 / 13 , 

aai = x (m^uj- 3 )^ 2 ^ 1 = xm 10 / 13 c^ 3 / 13 ,(30) 

r = ( rh W+i LU ™+^f +2fSrl = m 5 / 13 LU 18 / 13 . (31) 

Assumption of an optically thin gravitoturbulent disk is 
valid provided that a condition opposite to ([27)) is satis- 
fied. In particular, r < 1 for all CI < fit (or oj < 1) only 

if M < M f (or m < 1). 
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Properties of optically thin gravitoturbulent disks ex- 
hibit more moderate variation with r than in the opti- 
cally thick case. Indeed, when r < 1 and f3 — 2 disk 
temperature and aci vary with distance quite slowly, 
T oc r~ 3 / 13 and olqi oc r 9 / 26 . As a result, optically thin 
gravitoturbulent regime can be realized in a rather ex- 
tended region of the disk. 

As mentioned after equation (fT3"|) expressions (f2"3"|) - ([2"6")) 
and (j28|) - (f3~lj) also apply to disks in which M varies with 
distance. In this situation m = to(cj) and the dependence 
of various disk properties on uj can be directly obtained 
by plugging the expression for m(w) into these equations. 

3.3. External irradiation. 



According to equations (|24[) and (j2"9"|) temperature of a 
gravitoturbulent disk decreases with r. At some point ex- 
ternal irradiation becomes more important for the ther- 
mal balance than the internal gravitoturbulent dissipa- 
tion. Here we want to investigate a transition from a self- 
luminous disk (heated only by internal dissipation) to an 
irradiated regime. We assume that disk is illuminated by 
external radiation at temperature To. Given that most of 
the disk material is concentrated in its midplane region 
the gravitational stability of the disk is going to be deter- 
mined by the temperature T near the midplane and the 
condition for marginal stability is still given by Q rs Qo 
with Q defined in equation {1} . Below we consider sepa- 
rately the optically thick and optically thin cases. 

3.3.1. Optically thick case. 

Irradiated optically thick disk appears as an extension 
of a self-luminous optically thick disk. A specific location 
at which such transition occurs can be found by equat- 
ing T given by equation (|24[) to Tq. Angular frequency 
corresponding to this transition is 



LOT 



-1/3 [ ^0 

Tf 



(7-2/3)/6 



(32) 



By construction r at this location must be greater than 
unity, which with the aid of equation (|26|) yields the fol- 
lowing constraint on the optically thick transition to an 
irradiated disk: 



±£ > ^(5+2/3)- 

Tf 



(33) 



If r > 1 then the midplane temperature of the disk is 
given by the solution of the vertical radiative transfer in 
diffusion approximation as 



T w T n 4 



Mil 2 



rjT- 



(34) 



where 77 is a factor of order unity which absorbs our igno- 
rance of how the dissipation rate is distributed across the 
vertical thickness of the disk. Also, from energy conser- 
vation the temperature at the disk surface (photospheric 
temperature T p h) is 



1 ph — 1 Q 



um 2 

87TCT 
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These two relation allow us to distinguish three levels of 
the importance of irradiation. 



First, when 



ctT 4 < MQ 2 



(36) 



irradiation is so weak that it does not play any signifi- 
cant role even in the photosphere of the disk. This cor- 
responds to the case of a self-luminous disk which was 
covered in Sj3j 
Second, when 



MS} 2 < crT n 4 < tMVI 



(37) 



external irradiation keeps disk surface temperature at 
the level of To and creates a roughly isothermal gas layer 
underneath the surface. Nevertheless, the midplane tem- 
perature is still set predominantly by the internal dissi- 
pation and the outgoing flux ~ MQ 2 is still given by 
F = oT^/t, according to flS} and (37]). For that rea- 
son this irradiation regime also corresponds to the case 
of self-luminous disks explored in |J3] with some minor 
corrections having to do with the fact that the surface 
disk layers are hotter than they would have been in the 
absence of irradiation. 
Finally, when 



rMn 2 < o-t£ 



(38) 



irradiation is so strong that even the midplane tempera- 
ture reaches To and the disk becomes vertically isother- 
mal. This regime is different from the case of a self- 
luminous disk considered before since disk can no longer 
regulate its thermal state - its midplane temperature is 
fixed at T sa Tq. Assuming that disk is gravitoturbulent 
under strong irradiation we find from equation (jSJ) that 
its surface density must scale linearly with f2: 



1/2 



" «G Qe - - (f, ) »■ < 39 > 

where c Si o is the sound speed corresponding to tempera- 
ture To. According to this result E is independent of M 
in irradiated gravitoturbulent regions. 

Under strong irradiation one can no longer use the ex- 
pression ([4]) to derive aci- However, we know E from 
(|3"9"|) . so from the steady-state condition M = 37r^E and 
definition v = ac^/fl one finds that 



OiGI 



Qo GM 



-3/2 



(40) 



This result implies that ag; is constant in the strongly 
irradiated part of the disk. Note that analogous result 
has been previously obtained in Goodman (2003). 
We can also calculate the run of the optical depth 



T= CaKoT °n 

ttGQ 



Tf) 



0+1/2 



(41) 



which shows that r decreases as r~ 3 / 2 . This result com- 
bined with equation ([34]) also demonstrates that as r 
increases the relative contribution of the internal dissi- 
pation to the midplane temperature rapidly goes down 
since it is proportional to tVI 2 oc ft 3 . 

3.3.2. Optically thin case. 
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Optically thin irradiated disk can emerge as an exten- 
sion of an optically thin self-luminous disk, which hap- 
pens when T given by equation (|29[) becomes equal to 
Tq, or at 



LOT 



(9+2/3)/2 



(42) 



Such an optically thin transition is possible when an in- 
equality opposite to ([33]) is satisfied. 

An irradiated optically thin disk can also appear as a 
continuation of the irradiated optically thick disk con- 
sidered in §3 3.21 According to equation (|4Tj) this t = 1 
transition occurs at 
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(43) 



Optically thin disk is roughly isothermal vertically and 
its thermal balance requires 



T 



MCr 



V- 



(44) 



which is different from equation (|34|) by a factor r _1 
in the second term on the right-hand side accounting 
for the inefficiency of radiative cooling and absorption in 
the optically thin disk. There are two obvious regimes 
to consider. First, when 



<tT 4 < MQ 2 /t, 



(45) 



irradiation does not affect disk properties and we go back 
to the case studied in i )3 3.21 Second, when the condition 
opposite to ([4"5]) is satisfied irradiation sets the disk tem- 
perature. In this case all results (except for the equation 
[44] different from [34]) obtained in the optically thick 
irradiated case - equations ([39]) -([41 ]) - remain valid since 
in deriving them we did not use any assumptions about 
the value of r. 

The most important result regarding externally irradi- 
ated constant M disks is that they can remain gravitotur- 
bulcnt independent of their optical depth and that their 
effective viscosity aai remains constant. If the back- 
ground viscosity does not dominate angular momentum 
transport at least in some self-luminous parts of the grav- 
itoturbulent disk it will not dominate the transport in the 
irradiated part either since (Xgi is constant there. Thus, 
torque needed for transporting mass through the disk 
must be due to the GI. 

3.4. Fragmentation limit. 

When t coo i becomes comparable to fl^ 1 disk can no 
longer sustain the gravitoturbulence and has to fragment 
into bound objects (Gammie 2001). As mentioned before 
fragmentation condition Qt coo i < 1 can be recast in terms 
of the aai threshold according to equation (TH)]) . This 
formulation now allows us to directly apply our results 
for olgi derived in previous sections. 

In particular, self-luminous optically thick gravitotur- 
bulcnt disk starts fragmenting at 



LO < OJfrag 



,(4-2/3)/9 



(46) 



which follows from demanding oiqi given by equation 
([23]) to be larger than \ - the critical value of a needed 
for fragmentation. Radius r/ rag at which fragmentation 



first occurs is given by r/ rag = r $rn~^' 1 ~^l T! and is 
rather insensitive to either (3 or tin, so that fragmentation 
always occur not too far from r/. It is rather interest- 
ing that for /3 = 2, corresponding to the low tempera- 
ture dust opacity the location of fragmentation edge in 
the optically thick limit is completely independent of m: 
fragmentation occurs exactly at £1 = Qf. This fact has 
been first noticed by Matzner & Levin (2005). Clearly, 

constant M self-luminous gravitoturbulent disk cannot 
be fed by a source located outside v^ ra g. 

According to equation (|2T[) ujf rag corresponds to the 
optically thick part of the disk only if m > 1. Thus, 
whenever m > 1 the gravitoturbulent disk stays optically 
thick all the way to the fragmentation edge located inside 
of Tf. Then external disk feeding must necessarily occur 
interior to rj. 

In the optically thin case one finds from equation (|30p 
that 

co frag = rh^y\ (47) 

Fragmentation boundary lies in the optically thin part 
of the disk only if m < 1, in which case it is located 
outside of the fiducial radius r/. Fragmentation radius 
r frag = r/m~( 2/3+6 )/ 9 is a rather sensitive function of 
m: for (3 = 2 and M = 10~ 7 M Q yr" 1 « 0.014 M f 
(rather typical value of M in protoplanctary disks) one 
finds rf r ag ~ 10 2 rf. Thus, in the optically thin regime 
fragmentation can be pushed out to large distances by 
reducing m (for the just used values of j3 and M and 
M* = Mq one finds rf rag w 10 4 AU), but it still cannot 
be avoided if the disk is self-luminous. 

Everything we said before regarding fragmentation ap- 
plied to self-luminous disks. If the disk is stable against 
fragmentation all the way to the point where its tem- 
perature is determined by external irradiation then in 
the irradiated region aai is constant and given by equa- 
tion (|40p . Since oiqi must attain this value somewhere 
near the transition to irradiated regime and the disk is 
assumed to be non-fragmenting there (i.e. aai ^ X a t 
w ~ wt) we may conclude from equation ([40]) that the 
disk is going to remain in a gravitoturbulent state stable 
against fragmentation as long as it is externally irradi- 
ated and 

3/2 



< 



X 



(48) 



This is a rather interesting conclusion since it implies 
that a sufficiently low M disk can in principle be stably 
fed by a source of mass located at infinity. In partic- 
ular, for To « 10 2 K, when the low-temperature opac- 
ity still applies, one finds using estimates (|21[) that 
M < 10~ 4 M© yr -1 satisfies condition ([48]) in which we 
set x ~ 1 f° r simplicity. Thus, disks obeying the con- 
dition (|48p can transfer mass at a constant rate from 
very large distances despite being gravitationally unsta- 
ble, unlike the high M disks around quasars. Note that 
the criterion ([48]) is independent of either the opacity be- 
havior or the optical depth of the disk. A qualitatively 
similar conclusion about the stabilizing role of irradia- 
tion has been reached in Matzner & Levin (2005) and 
Cai et al. (2008). 

The argument based on equation (|48[) may not be with- 
out caveats. The value of aai given by equation ([40]) is 
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derived based on the cooling rate equal to the energy 
production rate ~ Mil 2 due to the dissipation of grav- 
itoturbulence. This rate is much lower than the (ex- 
ternal heating) rate ctTq x min(l,r) at which the disk 
would cool if irradiation were suddenly switched off or 
if the GI in irradiated disks were capable of producing 
surface density perturbations of order unity. One can 
easily show that Qt coo i based on such a fast cooling rate 
does not stay constant in irradiated part of the disk but 
steadily increases. The tricky question is the following: 
which cooling rate should determine the ability of the 
disk to fragment? This issue can be settled satisfactorily 
via careful numerical simulations of strongly irradiated 
gravitoturbulcnt disks, something that has not yet been 
done. But given that in accretion disks which are close to 
the steady state angular momentum transfer (determin- 
ing M) must be uniquely related to the energy dissipation 
rate we feel that it is more likely for the fragmentation 
condition to be determined by the criterion (|48[) rather 
than by the much shorter cooling time set by the irradi- 
ation heating rate. Our subsequent consideration will be 
based on this assumption. 

3.5. Background viscosity. 

If the angular momentum transport caused by gravi- 
tational torques becomes weak some other mechanisms 
may start providing effective viscosity. Here we assume 
that in addition to gravitational torques disk also pos- 
sesses some background viscosity a v due to e.g. the 
magneto-rotational instability (MRI). Normally one ex- 
pects a v <C 1 so that this background viscosity would 
become significant only when aci gets rather small. 

If, as expected for the dust opacity, a is between 
and 2 then both in the optically thick and optically thin 
regimes aai decreases as lo increases. In other words, 
gravitoturbulent disk becomes less "viscous" as the dis- 
tance to the center decreases. With this in mind we 
find that in the optically thick case the background vis- 
cosity would dominate over the gravitoturbulent torque 
(i.e. a v > o>gi) at 

w>w„« (af- 7 m 4 - 2 ^) 1/9 , r>l. (49) 

According to equation (|26|) this critical angular fre- 
quency corresponds to optically thick regime only if 

rh > a (10+4«/(7+4«_ 

In the optically thin case background viscosity regu- 
lates disk at 

w>w„« (a-^- 9 m^ +6 ) 1/3 , r<l, (50) 
and this io v corresponds to r < 1 provided that rh < 

(10+4/3)/ (7+4/3) 

Inside the region where the background viscosity dom- 
inates (at lo > u>„) equations governing disk structure 
change. Previously, when considering the gravitotur- 
bulcnt transport, we did not have a constraint on a 
but instead had a relationship between T and X in the 
form of equation (|9|), arising from the requirement of 
the marginal gravitational instability. However, with the 
background viscosity dominating the angular momentum 
transport a is constrained to be equal to a„, which leads 
to an overdetermined system of equations if we also try 
to keep the condition Q = Qq. This contradiction is 



naturally avoided by dropping the latter constraint, i.e. 
allowing the disk not to be marginally gravitationally 
unstable when a — a v . 

In Appendix [X] we present E, T, and r behaviors in 
accretion disk with the dominant background viscosity 
a v and opacity in the form ( [6]). If we now use these 
expressions [equations (|A1|) . (|A2[) , (|A5[) . and (|A6|) ] to 
calculate Toomre Q in the self-luminous viscous part of 
the disk at u> > lo v we find that 

Q/Q = (^) 9Al °- 2 «, t>1 (51) 
in the optically thick case and 

Q/Qo = (^M,) 3/(6+2/3) , t<1 (52) 

in the optically thin case. Apparently, for any reason- 
able dust opacity behavior Q starts to deviate from its 
marginal stability value Qo towards higher values at the 
transition from the gravitoturbulent to viscous regime 
(at lo = io v ). Thus, accretion disk cannot be gravitation- 
ally unstable if its angular momentum transport is not 
dominated by gravitational torques but is rather due to 
some other form of effective viscosity. This consideration 
demonstrates explicitly how the transition between the 
gravitoturbulcnt and viscous parts of the disk occur. 

Note that aai given by equation (|4"0")) is independent 
of lo. This implies that if 

/ rp \ 3/2 

fh < a v \ jrj , (53) 

accretion disk remains viscous everywhere and does not 
transition into the gravitoturbulent state at all. Thus, at 
low enough m a gravitationally stable viscous disk can 
transport matter all the way to the central object from 
infinitely large distances. Condition (|53|) is more restric- 
tive than the constraint (|48|) which delineates regime in 
which the gravitoturbulcnt constant M disk can transfer 
material from infinitely large distances without fragmen- 
tation. 

4. DISCUSSION. 

Our results derived in ij3] apply to a variety of situa- 
tions in which gravitoturbulent disks can exist: they can 
be optically thick or thin, fragmenting or not. Here we 
classify different states in which accretion disks can be 
found according to the values of M, background viscos- 
ity a v and irradiation temperature To. We also describe 
applications of these results to real astrophysical systems 
and discuss their connection with the work of others. 

4.1. Separation of different regimes. 

As we saw in previous sections different parts of ac- 
cretion disks can be characterized by different states — 
gravitoturbulcnt, viscous, irradiated, etc. Transitions be- 
tween these regimes depend on M, To, and a„. Our re- 
sults allow us to single out three important cases for the 
different transition topologies. 

When external irradiation is strong, 

To/Tf > 1, (54) 

the phase space of possible regimes can be represented by 
Figure [TJ Given that in the case of opacity dominated 
by cold dust T/ is just slightly higher than 10 K (see 
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Fig. 1. — Plot of different possible states in which an accre- 
tion disk can be found as a function of mass accretion rate M 
and angular frequency Q (corresponding dimensionless quantities 
fn = M/Mf and 10 = fi/f2/ are shown on the upper and right 
axes). Disk has a non-zero background viscosity «„ = 0.003 and 
is externally irradiated at a temperature To = 37/ £S 35 K, cor- 
responding to condition II 541 1. Shading indicates viscous regions 
(slant solid shading), irradiated regions (horizontal dotted shad- 
ing), and region where the disk must fragment (vertical solid shad- 
ing). Gravitoturbulent, self-luminous region is unshaded. Different 
curves separate regions with distinct physical conditions and are 
marked on the plot (see text for more details). Coordinates of the 
points where th ese curves cross in the (rh,u) coordinates are given 
by expressions (IB 1 ft . 



cq. [5T]) it is clear that the condition lj5"i)) should apply 
to virtually all accretion disks in the Universe as even 
the lowest measured temperatures encountered in some 
dense molecular cores are not very different from 10 K 
(Di Francesco et al. 2007). 

Using Figure [T] one can examine different regimes (in- 
dicated by shading) which are relevant for a given M. 
Curves 4 separating different regimes correspond to var- 
ious critical transitions: a = a v (dotted lines) im- 
plies a transition from viscous (above and to the left of 
this curve, slant solid shading) to gravitoturbulent state, 
a = 1 (solid lines) describes the onset of fragmentation 
(below and to the right of this curve, vertical solid shad- 
ing), T = T is a (long-dashed) line below which external 
irradiation starts to dominate disk structure (region with 
horizontal dotted shading), r = 1 curve (short-dashed) 
shows a transition from an optically thick (above this 
curve) to an optically thin disk. Coordinates (in phase 
space rn,uj, upper and right axes) of the points A, B, 
C, D where these curves cross are given by expressions 
(jBip in Appendix [Bl Unshaded region corresponds to a 
self-luminous, gravitoturbulent disk. A constant M disk 
corresponds to a straight vertical line in Figure [T] cutting 
through different regimes. 



4 Equations for these curves can be found in £j i|3 3.H3 3.51 and 
Appendix lAl 
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Fig. 2.— Same as Figure[T]but for a v = 0.003 and 7n = 0.77/ « 
8 K, corresponding to the condition (1551 ) . Coordinates of t he p oints 
where different curves cross are given by the expressions |)B2|I . 



The numerical values on the left and lower axes of this 
Figure correspond to a particular choice of a v — 0.003, 
To/T f = 3 (T = 35 K) and opacity in the form 0. 
As described in %T3~5l at the very low M < 10 -7 M Q 
yr _1 disk is always viscous and gravitationally stable, 
even at arbitrarily large distances from the central ob- 
ject. Above this value of M the disk must be gravito- 
turbulent within some range of distances. For 10 -7 M Q 
y^ 1 1$ M < 4 x 10" 5 M yr -1 gravitoturbulent region 
extends from ~ 10 -10 — 10~ 9 s _1 all the way to very 
large distances (where ft — > 0). However, for M above 
M > 4 x 10~ 5 M Q yr" 1 gravitoturbulent disk fragments 
at O « 1.4 x 10~ 10 s _1 (in the optically thick case with 
opacity characterized by f3 = 2 fragmentation occurs ex- 
actly at VI f, see H3 3.4|) so that a constant M disk can be 
maintained only interior to this point. For any M > 10 
Mq yr^ 1 the gravitoturbulent angular momentum trans- 
port becomes weak at large enough SI so that the back- 
ground viscosity starts to determine disk properties at 
small radii (upper part of the plot, at Jl > 10~ 9 s _1 in 
this particular case). 

As an example, let us use Figure [T] to figure out where 
the transitions between different regimes occur in a disk 
with M = 1Q- 5 M Q yr" 1 , a v = 0.003, T /T f = 3 and 
opacity in the form (J7|. Atf2>10 -9 s~ 1 disk is grav- 
itationally stable and self-luminous, angular momentum 
is transported through the disk by background viscosity. 
For 2 x 10~ 10 s _1 < fl < 10~ 9 s _1 disk is gravitoturbu- 
lent and self-luminous. For £1 < 2 x 10~ 10 s _1 extending 
to infinity disk is gravitoturbulent and externally irradi- 
ated. Disk is optically thick for all fl > 10~ n s _1 and 
optically thin for smaller Q. This example clearly shows 
the complexity of possible states in which a given accre- 
tion disk can be found at different distances from the 



central object. 

Topology of the phase space of possible disk states 
changes when the temperature of external radiation field 
becomes smaller than Tf . This situation is not very typi- 
cal (remember that the CMB temperature in the present 
day Universe which would set a temperature floor is not 
too far from Tf) for astrophysical objects although tem- 
perature at the level 6 — 7 K has been found in centers 
of some dense molecular cores well shielded from exter- 
nal starlight, so that situations in which To < Tf can 
be relevant for very young protoplanetary disks (Class 
objects without central heating source) forming inside 
very dense cores. 

For completeness we display the different topologies of 
disk states possible for To < Tf. There are two cases. 
First, for 



v 2/(7+4/3) 



< T Q /T f < 1 



(55) 



the separation of different regimes in the M, f2 space is 
represented by Figure [5J One can see from this Figure 
that the range of M in which there is a gravitoturbulent 
region extending to infinity shifts to lower M (Figure [5] 
is drawn for the same a v — 0.003 as Figure [1] but for 
T = 0.7T f w 8 K). 

Second, for very low irradiation temperature 

T /Tf < oftV+W (56) 

the separation of different regimes is shown in Figure [3J 
This plot again assumes a v = 0.003 while T = 0.3T/ « 
3.5 K. Such low temperature is hardly achievable in re- 
ality as it is just slightly higher than the current CMB 
temperature. Note that the topology of different regimes 
in Figure [2] is intermediate between the cases shown in 
Figures [1] and [3] — it looks like the former at low M and 
like the latter at high M . 

Despite some differences in details of how exactly the 
transition from viscous to gravitoturbulent, to irradiated 
states occurs at intermediate values of M, Figures Q][3] ex- 
hibit four invariant properties of accretion disks: (1) disk 
is dominated by background viscosity for all O at very 
low M [for m < a u (T /Tf) 3 / 2 , see Appendix IB] and non- 
zero To, (2) for all M disk is dominated by background 
viscosity at high enough il, (3) at intermediate values of 
M [for a u (T /T f ) 3 / 2 < rh < {T Q /T f f/ 2 ] disk possesses 
a gravitoturbulent, externally irradiated region that ex- 
tends to arbitrarily large distances, (4) at large values of 
M [for rh > (To/T/) 3 / 2 ] disk has a gravitoturbulent re- 
gion within a finite range of distances but must inevitably 
fragment at some large distance. Similar qualitative con- 
clusions hold also for opacity behaviors different from ((7]) . 

We should note here that calculations presented in this 
section rely on our use of viscosity ©-([7]) throughout the 
whole region of M, f2 phase space that we consider. In 
reality, at high M and disk temperature should exceed 
10 2 K at which point icy grains sublimate leaving metal- 
silicate grains as a source of opacity. This latter opacity 
source while still being in the form ([6]), is characterized by 
smaller values of k at the same temperature and (3 1/2. 
This is likely to quantitatively (but not qualitatively) 
affect results presented in Figures [T][3] at high M and il. 
To not complicate things further here we do not attempt 
to self-consistently describe transitions between different 
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Fig. 3.— Same as Figure [TJ but for a v = 0.003 and T = 
0.3Ty Ri 3.5 K, corresponding to the condition 1561 . Coordinates of 
the p oints where different curves cross are given by the expressions 
(tB3l . 
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opacity regimes but rather display a qualitative pict 
for a single opacity law. 

4.2. Applications. 

Our results can be applied to understanding the prop- 
erties of the outer, cold parts of realistic accretion disks. 
In particular, we address three important issues. 

First, a possibility of giant planet formation by GI in 
protoplanetary disks has been discussed since Cameron 
(1978). In this context it is interesting to ask under which 

conditions a constant M protoplanetary disk would be 
prone to fragmentation into gravitationally bound, self- 
gravitating objects. Our results described in t j44.ll can 
directly address this issue. Indeed, conditions used in 
producing Figure [IJ namely a = 0.003 and To w 35 K 
are quite typical for the outer regions of protoplanetary 
disks, beyond ~ 100 AU from the central star. One one 
hand, the disk surface density there is low enough (see 
below) for the cosmic ray ionization to stimulate MRI 
operation which gives rise to a v at the level of ~ 10~ 3 — 
10~ 2 . On the other hand, outer regions of protoplanetary 
disks are warmed up by radiation of either the parent 
star or the neighboring stars at the level of several tens 
of K. Thus, the situation represented in Figure [T] can be 
directly used for understanding the properties of external 
parts of protoplanetary disks. 

What is obvious from this Figure is that giant planet 
formation by gravitational instability in constant M 
disks can take place only beyond « 120 AU which is the 
distance from a 1 Mq star at which Q = ilf — remember 
that in the optically thick 5 gravitoturbulent disks frag- 
mentation occurs at this specific value of fl (Matzner & 

5 That the disk with the assumed values of a v and To is optically 
thick at fragmentation boundary is evident from Figure [T] 
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Levin 2005), see equation (146)) . Thus, planets produced 
by gravitational instability should be born far from their 
parent stars although one cannot exclude their subse- 
quent migration to shorter periods. 

Another obvious constraint on planet formation that 
follows from Figure [1] is that M must be pretty high 
at the location where disk fragments and planets form. 
Indeed, one can easily see that fragmentation is possible 
only if M locally exceeds 10~ 5 M© yr -1 . At lower M 
disk maintains itself in a gravitoturbulent state (or even 
being kept gravitationally stable by its own background 
viscosity at very low M) without fragmentation even very 
far from the star. 

Accretion rates in excess of 10 -5 M Q yr _1 are atypical 
for mature T Tauri disks (Gullbring et al. 1998). How- 
ever, they may have been present at the very earliest 
stages of star and disk formation when the material from 
collapsing protostellar envelope rains down onto the disk 
at a very high rate, possibly exceeding 10 -5 M Q yr _1 
in some locations. Such disks are likely not to have M 
constant through their whole extent but as we discussed 
in S}3] our results are still applicable 6 even in this more 
complicated case as long as M is specified locally. 

Second practical issue that we are going to address has 
to do with the feeding of supcrmassive black holes in 
centers of galaxies. It has been known for a long time 
that the outer parts of quasar disks must be gravita- 
tionally unstable which was always raising a question of 
how gas is transported to the black hole from large dis- 
tances. Our results demonstrate that as long as M is 
not very high the GI is not going to impede mass trans- 
fer through the disk since for low enough M, namely for 
M < M f (T /T f ) 3 / 2 , see equation (g5]), a gravitationally 
unstable disk can persist in a gravitoturbulent state at 
arbitrarily large distances from the central object. How 
high M can be carried through a gravitoturbulent disk 
globally thus depends only on the level of external irra- 
diation. 

Radiation fields in galactic nuclei due to circumnuclear 
stars are expected to be quite intense giving rise to To at 
the level of tens to hundreds of K. Assuming To = 100 K 
(which is a radiation field slightly more intense than that 
expected in the Galactic Center) we find that gravitotur- 
bulent disk can transport mass from infinity to the black 
hole as long as M < !0 3 M f « 10~ 2 M Q yr _1 . This is 
about 10% of the Eddington rate (for radiative efficiency 
of 10%) for the 4 x 10 6 M Q black hole in the center of our 
Galaxy, which is quite significant given that the Bondi 
accretion rate of this object is < 10~ 5 M© yr _1 (Baganoff 
et al. 2003). Thus, irradiated gravitoturbulent accretion 
disks provide a natural way of continuous feeding at least 
some (not too massive) central black holes at reasonable 
rates by gas transported from very large distances. 

However, it is also clear from our results and Figure 
Q] that accretion disks around more massive black holes 
(> 10 6 M ) consuming mass at rates close to Edding- 
ton cannot remain gravitoturbulent out to very large dis- 
tances — for any reasonable level of external irradiation 

6 A disk with M varying with distance would not correspond to 
a straight vertical line in Figures 11131 like a constant M disk would 
but must follow a more complicated path determined by a specific 
dependence of M on r (or, alternatively, f2). 



disk must fragment at some point, and the most distant 
possible location of the fragmentation boundary corre- 
sponds tof2 = r2/«1.4x 10~ 10 s _1 (for the case of dust 
opacity in the form [7]). Transport of gas from beyond 
this distance is still an open issue discussed by Goodman 
(2003). 

Finally, we briefly discuss the origin of young stellar 
disk around a supermassive black hole in the center of our 
Galaxy. Inner parsec of the Galaxy is known (Paumard 
et al. 2006; Lu et al. 2009) to contain at least one disk 
of young (6 ± 2 Myr old) massive stars spread between 
0.04 pc and 0.5 pc. To explain formation of these stars so 
close to the black hole where they are subject to action of 
its strong tidal field a fragmentation of a gravitationally 
unstable disk has been proposed (Levin & Beloborodov 
2003; Levin 2007). Such an event can in principle hap- 
pen both in a (quasi-)steady state disk like the one we 
considered in this work or in a short-lived promptly frag- 
menting (on a dynamical time scale) disk-like structure 
that may arise as a result of molecular cloud collision 
with the black hole (Wardle & Yusef-Zadeh 2008; Bon- 
nell & Rice 2008; Hobbs & Nayakshin 2008). Here we 
try to constrain the first possibility. We will assume that 
disk was illuminated by surrounding stars which kept To 
at the level of tens of K so that its opacity law was given 
by equation ([7]). Distance rt at which fragmentation 
would occur in an optically thick disk around the 4 x 10 6 
M Q black hole is w 0.1 pc which is within the span of 
the observed stellar disk. 

For a long-lived disk to start fragmenting a variabil- 
ity of some of its properties must be taking place. One 
possibility is an increase of M which can bring a gravito- 
turbulent irradiated disk extending out to large distances 
across a fragmentation threshold (see Figure[T|). But then 
the outer regions of the disk (at r ^> r/) where M has 
increased would immediately fragment and it is not at all 
obvious that mass would be transferred inward increas- 
ing M at small radii (at r ~ rf) where stellar disks are 
observed. In principle, disk could be not an accretion but 
a spreading (Pringle 1991) disk formed as a result of a 
dense molecular cloud disruption at very small distances. 
However, in this case an increase of M is unlikely while 
a decrease of M typical for spreading disks would only 
stabilize the disk against gravitational fragmentation. 

Another (probably less likely) possibility is a reduction 
of external irradiation which can cause fragmentation of 
even a constant in time M disk as soon as the condition 
(j48|) gets violated. The problem with this scenario is that 
it may then be difficult to explain the existence of stars at 
0.04 pc which is significantly closer to the black hole than 
the minimum radius Tf at which fragmentation occurs in 
a cold disk. In principle, Qf can be increased by lowering 
disk metallicity which affects kq. However, according to 
equations (fT7| & (|22|) moving 77 from 0.1 pc to 0.04 pc 
would require reducing kq by a factor of 60 compared to 
the value in |7|). implying extremely sub-solar metallicity 
in the disk. 

More generally, star formation in a nearly-Keplerian 
disk does not naturally explain rather significant eccen- 
tricities of disk stars (Bartko et al. 2008) and the possible 
presence of a second disk component. Thus, it seems un- 
likely (although not completely impossible) that stellar 
disks were formed by gravitational fragmentation of a 
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long-lived gaseous disk. Scenario of a prompt fragmen- 
tation of a tidally disrupted molecular cloud advanced 
by Wardle & Yusef-Zadeh (2008), Bonnell & Rice (2008) 
and Hobbs & Nayakshin (2008) presents a more attrac- 
tive possibility. 

4.3. Comparison with previous studies. 

The first investigation of self-gravitating accretion 
disks has been undertaken by Paczynski (1978a), later 
followed by Paczyncki (1978b) and Kozlowski et al. 
(1979). These early numerical modelling efforts concen- 
trated on studying hot quasar disks around the super- 
massive black holes with accretion rates close to the Ed- 
dington rate. Some of the important ingredients of these 
models have been the inclusion of the radiation pressure 
(neglected in our case) and use of high-T opacities, which 
makes comparison of these calculations to the results of 
our study rather difficult. 

Another investigation of a quasi- viscous evolution of a 
self-gravitating disk driven by gravitational torques has 
been done by Lin & Pringle (1987). In their work, based 
on rather general arguments, a specific model for the 
viscosity due to the disk self-gravity has been assumed, 
namely v oc £ 2 r 6 f2 (or a ~ Q~ 2 > 1). The correct pre- 
scription (|2|) is in general quite different from this naive 
anzatz precluding direct comparisons with our results. 
Also, the assumption Q ~ 1 has been relaxed in Lin 
& Pringle (1987) allowing Q to drop significantly below 
unity. 

More recently Goodman (2003) has analytically inves- 
tigated properties of Q ss 1 regions of constant M quasar 
disks, again taking into account radiation pressure. This 
study has a lot in common with our work with the major 
difference being that Goodman (2003) left a to be a free 
parameter while we self-consistently calculate its value 
using a prescription ^ . Some of the results derived in 
Goodman (2003) have been retrieved in our study. 

Subsequently, Rafikov (2005, 2007) has looked at the 
properties of gravitationally unstable disks which are ca- 
pable of forming giant planets by direct fragmentation. 
Such disks must simultaneously fulfill two constraints: 
Q = Qo and o>gi — X- These assumptions were used to 
set stringent constraints on the properties of disks that 
are able to form planets by GI. Clearly, such disks do 
not have constant M in general. Our current assump- 
tions are different in that we assume only Q = Qq and 
fix M at some value (not necessarily constant, see the 
discussion after equation |13| ) at every point in the disk, 
but aci is then calculated self-consistently. Matzner & 
Levin (2005) and Levin (2007) have also studied cold 
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Q = 1 disks with opacity dominated by dust specifically 
looking at the conditions necessary for fragmentation. 
We successfully reproduce some of their results such as 
the location of the fragmentation boundary in the opti- 
cally thick regime and the importance of irradiation for 
stabilizing the gravitoturbulent disk against fragmenta- 
tion at large distances. 

Finally, Terquem (2008) & Zhi et al. (2008) have con- 
structed global numerical models of protoplanetary disks 
accounting for the possibility of GI in some parts of the 
disk. These studies pay special attention to the presence 
of the so called "dead zones" (Gammie 1996) — disk 
regions where MRI cannot operate because of low ion- 
ization. Our analytical calculations do not account for 
the existence of such regions, nevertheless they provide 
good foundation for understanding numerical results of 
Terquem (2008) & Zhi et al. (2008) in the outermost re- 
gions of their disks. 

5. CONCLUSIONS. 

We have explored the properties of marginally gravi- 
tationally unstable accretion disks using a realistic pre- 
scription for the angular momentum transfer driven by 
the gravitational torques. We self-consistently derived 
scalings of important disk variables such as surface den- 
sity and temperature in both the optically thick and 
thin regimes. We also accounted for the possibility of 
disk having some background viscosity, e.g. due to MRI, 
and demonstrated that in this case a gravitoturbulent 
disk inevitably switches to an ordinary viscous disk at 
small radii. Another important ingredient of our study 
is the inclusion of possible external irradiation of the disk 
(e.g. by central object or nearby stars). We have demon- 
strated that for low enough M external irradiation helps 
to stabilize the disk against fragmentation all the way to 
infinity providing means of mass transport to the cen- 
tral object from very large distances. At extremely low 
M disks have been shown to never become gravitation- 
ally unstable because of background viscosity. At high 
M (the exact threshold depends on the irradiation tem- 
perature) fragmentation of the disk is inevitable at large 
distances. Results of this work apply to our understand- 
ing of the possibility of giant planet formation by GI and 
star formation in the Galactic Center, and to the problem 
of feeding the quasars. 
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APPENDIX 
PROPERTIES OF VISCOUS ACCRETION DISKS. 

When the angular momentum transport in the disk is dominated by some internal source of viscosity characterized 
by parameter a v one finds that in the optically thick case 

,3-4-1 (5-/9)-" 



Kq 



a 



4-/3 



3 KqM fj, 



(5-/3)" 



T 



(5-/3)- 1 



4-/3 



T>1, 



where we have dropped constant numerical factors of order unity. With these expressions one finds 

(5-/3)" 1 



meaning that viscous disk is optically thick only when 



w > 



1/3+3 



(2+2/3)- 



Whcn the condition (|A4[) is violated disk becomes optically thin and we find 



n 2+/3jy-3+/3 
"4+73 



k B 



4+/3' 



(3+/3)" 



2+/3 \ (3+/3)- 1 



4+/3 



(3+/3)" 



<JKq fi 
,1+0 x 2/(3+0) 



rr I \(3+/3) _1 



T < 1. 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 

(A6) 
(A7) 



Derivation of equations (|A1|) - () A7|) assumes that it is the internal viscous dissipation in the disk that sets its midplane 
temperature. However, analogous to the case studied in ^3 3. 31 one can consider a possibility that the disk temperature 
is set by external irradiation at the level of T = Tq even when its angular momentum transport has non-gravitational 
nature. From equation (|A2|) we find that an optically thick self-luminous viscous disk changes to an externally 
irradiated disk at a radius where 

(5-3)/3 



\m 2 ) \T f 



(A8) 



Analogously, equation (|A6[) predicts that an optically thin transition from a self-luminous to an externally irradiated 
disk occurs at the point where 

(3+13) 



In those regions where the disk temperature is fixed at the constant level To one can easily show that 

rilLO Tq 



OL v If 

a v \T f J 



(/3-1) 



Q/Qo 



in 



Tf) 



3/2 



(A9) 

(A10) 
(All) 
(A12) 
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Last equation implies that Toomre Q is constant in the irradiated viscous part of the constant M disk. That Q > Qq 
there can be easily seen from equation (140)1 which allows us to rewrite equation (|A12[) as Q/Qo = a^/aci and this 
ratio is > 1 since for a u to dominate over the gravitoturbulcnt torques a v > aci must be fulfilled. It also follows from 
equation (|A11[) that the r = 1 transition, if it occurs in the externally irradiated region of viscous disk, takes place at 

— £®H- (A13) 

PHASE SPACE PLOTS. 

The (rh,u) coordinates of various critical points in Figure [T] corresponding to Tq/Tj > 1 are 

/ / T \3/2 / T \(2-/3)/3\ / /t n3/2 , „ x -(/3+l/2)N 



TfJ v \T s j j 1 V \T f J '\T f 

rp \ 3/2 /rp x -(/3+l/2)\ f / rp \ 3/2 /r^ \ (2-,8)/3N 



C "U57J A^J j' l^ 1 (B1) 

The same for a^ /(7+4/3) < T /T f < 1 displayed in Figured 



A = \ a v — , a„ ' J — , B = \ a v — , — 



2/(7+4/3) 



and for Tq/T/ < a„ shown in Figure [3] 
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